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Abstract 

We study the manganese oxides from the viewpoint of the strongly corre- 
lated doped Mott insulator. The magnetic ordering and the charge trans- 
port are governed by the orbital degrees of freedom, and their dimension- 
ality is controlled by the anisotropic transfer integrals between the e g or- 
bitals. As x increases the magnetic structure is predicted to change as 
A ^ F — > A ^ C ^ G { F: ferromagnet, A: layered antiferromagnet, 
C: rod-type antiferromagnet, G: usual antiferromagnet), in agreement with 
experiments. Especially the orbital is aligned as d x i_ y 2 in the metallic A state, 
which explains the quasi 2D transport and no canting of the spin observed 
experimentally. 

Next we discuss the ferromagnetic state without the orbital ordering due 
to the quantum fluctuation. Here the interplay between the electron repulsion 
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U and the Jahn- Teller electron-phonon interation E'lr is studied with a large 
d model. In addition to this strong correlation, we propose that the dynam- 
ical phase separation could explain the specific heat as well as the various 
anomalous physical properties, e.g., resistivity, photo-emission, etc. 

I. INTRODUCTION 

The interplay between the magnetism and transport properties is one of the main issues 
in the physics of the strongly correlated electronic systems. Doped manganites Ri-^A^MnOs 
(R=La, Pr, Nd, Sm ; A= Ca, Sr, Ba) have recently attracted considerable interests from this 
viewpoint in connection to the colossal magnetoresistance (CMR). It has been believed that 
the most fundamental mechanism for this CMR is the double exchange interaction, which 
has been studied long time ago [0-01 • However recent systematic experimental studies have 
revealed the vital role of the orbital degrees of freedom, which is the subject of this paper. 

In the parent compound RMnOs, the valence of the Mn ion is Mn 3+ with the configura- 
tion d A . Due to the crystal field the five d-orbitals split into two e 9 -orbitals (x 2 — y 2 , 3z 2 — r 2 ) 
and three t 2s -orbitals (xy, yz, zx). According to Hund's rule the four electrons occupy three 
/f^g-orbitals and one of the e s -orbitals with parallel spins. Then the e 9 -electron has the or- 
bital degrees of freedom in addition to the spin ones. When the holes with concentration x 
are doped, the vacancies in the e 5 -orbitals are introduced which induces the ferromagnetism 
due to the double exchange interaction and the charge transport. However there are several 
experimtental facts which cannot be explained in the simple double exchange models. One 
is the rich phase diagrams which is studied theoretically in Section II and the other is the 
anomalous features in the ferromagnetic metallic phase described in Section III. All these 
issues are colsely related to the orbital degrees of freedom, which control the dimensioanlity 
of the electronic systems and whose quantum fluctuation gives rise to the orbital liquid and 
orbital Kondo effect. Section IV is devoted to the conclusions where the possibility and 
implications of the dynamical phase separation are mentioned. 
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II. PHASE DIAGRAM 



We set up the three-dimensional cubic lattice consisting of the manganese ions. Two 
kinds of the e g orbital( 7, 7') are introduced in each site, and the t 2g electrons are treated 
as the localized spin with S = 3/2. The Hamiltonian without the JT coupling is given as 
follows §!, 

H = E (tl/4 7a d jYrT + h.c.) 

<ijr'>,<r,7,7' 

+ UY, n i^ n Hi + U'J2 H ia n ib + 1 4aa4ba' d iaa' d iba 
i-y i i,a,a' 

+ j H j2s t -s^ + j s £ st 2s -s^ . (i) 

i <ij> 

is the operator which creates an electron with spin cr(=|, j) in orbital 7(= a, b) at site i, 
and Si is the spin operator for the e g electron defined by Si = \ J2aa'y dia-y^^'dia'-y- The first 
term describes the electron transfer between 7 orbital at site % and 7' orbital at the nearest- 
neighbor site j. The value of t]J is estimated by considering the oxygen 2p orbitals between 
the nearest Mn-Mn pair, and is represented by c]J to, where c]J is the numerical factor 
depending on the orbitals 0. t is estimated to be 0.72eV, which we choose the unit of energy 
below (to = 1) UJ]. The second line shows the electron-electron interaction terms where U, U' 
and / is the intra-, inter-orbital Coulomb interactions, and inter-orbital exchange interaction, 
respectively. This interaction can be rewritten as -a Et^iW + ^^"(Y)) — /3SiTj(r) 2 
0. Here the spin operator Si and the iso-spin operator Tj = | I] 77 / CT <iJ - 7 o ! 77 'C?i - 7 ' for the 
orbital degrees of freedom are introduced, and the two positive coefficients a and (3, which 
are defined by a = 2U/3 + U'/3 — 1/6 and ft = U' — 1/2 , represent the interaction to induce 
the spin and iso-spin moments, respectively. The last line is the sum of the Hund coupling 
and the AF interaction between the nearest neighboring t 2 g spins. Here we adopt the mean 
field approximation by introducing the order parameters (Si), (>5* 29 ), and (T/). These order 
parameters are determined to optimize the mean field energy at zero temperature. For both 
spin and orbital, the four types of the ordering are considered, that is, the ferromagnetic (F- 
type) ordering, where the order parameters are uniform, and the three AF-like orderings, i.e., 
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the layer- type (A- type), the rod- type (C-type) and the NaCl-type (G-type) AF orderings. 
Hereafter, types of the orderings are termed as, for example, (spin:C), and so on. 

In Fig. 1, we present the calculated phase diagram in the case of a = 8.1, (3 = 2.5, which 
is relevant to the actual manganese oxides. In each region, the spin structure is specified by 
the letter F,A,C,G and the orbital configuration by picture. The important results are the 
followings. 

[1] There occurs the crossover from the superexchange dominated region for small x (x < 0.3) 
to that of the double exchange interaction for large x (x > 0.3). The former region could 
not be reproduced by the analysis without Coulombic repulsion. In this respect, the spin F 
state up to x ~ 0.4 in the phase diagram is mostly due to the superexchange interaction. 
[2] Because of the anisotropy of the transfer integral between the e 9 -orbitals, the most natural 
orbital configuration expected from the double exchange interaction is (x 2 — y 2 ) with the 
two dimensional electronic dispersion. This corresponds to the A phase for 0.15 < x < 0.45. 
In this case the interlayer transfer is forbidden, and the spin A-type sturucture is stabilized 
by J s . No spin canting is expected in this case. 

[3] Similarly the spin C structure is accompanied by the 3z 2 — r 2 orbital which gives the 
ID-like dispersion along the c-axis. Hence the spin and orbital structure is closely correlated 
through the dimensionality of the electronic dispersion. From this viewpoint, the orbitals 
have more degrees of freedom in the ferromagnetic state. This suggests the large orbital 
fluctuation in the spin F state as discussed in the next section. 

We now discuss the comparison between the mean field phase diagram in Fig. 1 and the 
experiments. At the moment, a value of Js cannot be estimated accurately, but there are 
two rough estimates. One is from the Neel temperature T/v = 130/^ for CaMnOs (x = 1.0) 
||, which suggest Js = T/v/7.5 = 1.7meV = 0.0023to in the mean field approximation. The 
fluctuations lower Tjy, and hence increase the estimate for Js- Another estimate is obtained 
from the numerical calculations for LaMn03 (x = 0.0), which suggests Js = 8meV = O.Ollto 
||]. Although Js might depend on x in real materials, we tentatively fix Js to be 0.009to 
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represented by the broken line in Fig. 1. Then the spin structure changes as A — > F 



C — > G, as x increases, which is in good agreement with the experiments [1C-12 



Another implication to the experimental results is about spin:A phase appearing around 
x > 0.5. In Ndi_ x Sr x MnO"3, the ferromagnetic metallic phase is realized up to about x = 0.48 
and the CE-type AF structure with the charge ordering tunes up P,|T0"|. With further 
increasing of x, the metallic state with spin:A again appears at about x = 0.53, and the 
large anisotropy in the electrical resistivity is observed in this phase. The similar metallic 
phase accompanied with spin: A is also reported in Pri_ x Sr x MnO"3 jfl]]. The spin canting is 
not observed experimentally [pHTTfl as predicted above in this Spin:A metal, in contrast to 
the Spin: A in the small x region which has been discussed by de Gennes 0. 

Summarizing this section we have studied the phase diagram of Lai-^Sr^MnOa in the 
plane of x ( hole concentration ) and J s ( AF exchange interaction between the t2 9 spins) in 
the mean field approximation. The global features can be understood in terms of the inter- 
play between the superexchange interaction and the double exchange one which is consider- 
ably modified with taking the orbital degrees of freedom into account. The dimensionality 
of the energy band attributed to the orbital structure plays an essential role to determine 
the phase boundary. The orbital structure in the ferromagnetic phase is sensitive to changes 
of the carrier concentration and the energy parameters. This is because the spin structure 
does not forbid the kinetic energy gain for any direction, and the orbitals have more free- 
dom. This enhances the fluctuation of the orbital degrees of freedom. Actually the orbital 
fluctuation has been studied in the slave-fermion formalism, and it has been found that the 
orbitals remain a liquid down to the low temperature due to the fluctuation among the three 
directions of the x 2 — y 2 -type ordering || . What quenches the entropy of the system is the 
quantum tunneling between the different orbital configurations, which could be regarded as 
the orbital-Kondo effect. This will be studied in the next section. 
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III. ORBITAL FLUCTUATION 



Experimentally there are several anomalous features in the low temperature ferromag- 
netic phase below T c . 

[a] In the neutron scattering experiment no temperature dependent phonon modes have 
been observed The recent Raman scattering experiment also shows that the Jahn- 

Teller phonons ( especially their frequencies ) are temperature independent and insensitive 
to the ferromagnetic transition at T c [jnfl . [b] The photo-emission spectra show a small 
discontinuity at the Fermi edge even at T << T c , which suggests some interactions still 
remain strong there |15| . [c] The optical conductivity cr{uj) at T « T c is composed of two 



components, i.e., the narrow Drude peak (a; < 0.02eV) and the broad incoherent component 



extending up to u ~ leV [T6| . The Drude weight is very small, which seems to be consistent 



with the photo-emission spectra. 

[d] The low temperature resistivity p(T) can be fitted by 

p(T)=p + AT 2 , (2) 

where A is a large constant of the order of 500yuficm/K 2 ||17|| , again suggesting the strong 
electron correlation. 

[e] Contrary to the case of resistivity [d] , the coefficient of T-linear specific heat is very small 
with 7 = 2mJ/K, which violates the Kadowaki- Woods law for these compounds flUf] . 



Although [e] is difficult to reconcile with [a]-[d], we consider the latter as the evidence 
for the strong correlation and orbital fluctuation even at T T c . Because the spins are 
perfectly aligned at T <^ T c , the only remaining degrees of freedom are the orbital ones. In 



this section, we study a large d model |19[ for the ferromagnetic state including both the 



electron-electron interaction U and the Jahn- Teller coupling g ||20|| . This model takes care 
of the following two respects, (i) including the electron-electron interation, (ii) including 
the quantum fluctuations. Especially the latter is essential to describe the low temperature 



Fermi liquid state, which is described by the Kondo peak in the large d limit |19| . The strong 
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electron-electron interaction with the reasonable magnitude of the Jahn- Teller coupling ex- 
plains both the large isotope effect and [a]. Moreover, the features of strong correlation 
[b]-[d] are at least consistent with the large U picture although [e] still requires another 
physics, which we will discuss in the next section. We start with the Hubbard-Holstein 
model for the ferromagnetic state. 

h = - E C'4cia' + £/Ewa 

i,j,a,a' i 

- 9 E Qifa - raj) + ~ E(^ + MQ 2 Q 2 t ), (3) 

i i 

where the spin indices a ="[*, [ correspond to the orbital degrees of freedom as f = d x 2_ y 2 and 
1= dz z 2_ r 2, and the real spins are assumed to be perfectly aligned. We consider only one 
Jahn- Teller displacement mode Qi for each site, while there are two modes Q2iiQzi in the 
real perovskite structure . However this does not change the qualitative features obtained 



below. The Jahn- Teller mode is assumed to be an Einstein phonon with a frequency Q. The 
transfer integral t^' depends on a pair of orbitals (a, a') and the hopping direction. These 
dependences lead to the various low lying orbital configurations, and thus they suppress 
the orbital orderings in the ferromagnetic state. Actually there is no experimental evidence 
for the orbital orderings, e.g., the anisotropies of the lattice constants and/or transport 
properties in the ferromagnetic state. Then we assume that no orbital ordering occurs 
down to the zero temperature due to the quantum fluctuations, and the transfer integral 
is assumed to be diagonal for simplicity, i.e., tff = Uj5 aa '- This means that the ground 
state is a Fermi liquid with two degenerate bands at the Fermi energy. In order to describe 
this Fermi liquid state, we employ the large-c? approach where the model Eq.(3) is mapped 
to the impurity Anderson model with the self-consistent condition This approach has 



been applied to the manganese oxides to study Hund's coupling by Furukawa an d to 
study the additional Jahn- Teller coupling by Millis et al. |23|]. The action for the effective 
impurity model is given by 

rP rP 



S=f dr [ dr'cl{T)G,\T-T')c a {T') 
Jo Jo 



+ / dr 



Jo 



Un^T) ni (r) - #Q(t)(ti t (t) - n^r)) 



/3 , M 
ar — 
o 2 



(<9 t Q(t)) 2 + ^ 2 Q(t) 



, (4) 



where G 1 (r — r') is the dynamical Weiss field representing the influence from the sur- 
rounding sites. The self-consistency condition is that the on-site Green's function G(iu n ) = 
(G (i^ n ) _1 — S(ic<j n )) _1 should be equal to the Hilbert transform of the density of states 
D(e) as 

D(e) 



G{iuj n ) = de 
J 



(5) 



iuj n + ji — Tl(icu n ) — e 

We take the Lorentzian density of states D(e) = t/(it(e 2 + t 2 )) because there is no need to 
solve the self-consistency equation in this case, and the Weiss field is given by 



G 1 (iu n ) = iuj n + /! + it sign lu Ti 



(6) 



The only unknown quantity is the chemical potential /i, which is determined by the electron 
number. The determination of chemical potential requires some numerical calculations, but 
our discussion below is not sensitive to the location of the chemical potential. Then we take 
ji as the parameter of our model, and the problem is now completely reduced to that of a 
single impurity. Now we introduce a Stratonovich-Hubbard (SH) field £(r) to represent the 
Coulomb interaction U. Then the action is given by 

S = J2( iuj n + V + it sign u n )c ] a (u n )c a (u n ) 



dr 



A' 



£0 



M 



[(d T Q(r)) 2 + n 2 Q(r) 2 ] + C(t)Q(t) - (A£(r) + gQ(r)){n^r) - H {t)) 



where A = (t 2 + /i 2 )/t and £(r) is the test field to measure the phonon correlation function. 
At this stage the electron is coupled with the linear combination of the SH field £ and the 
phonon Q as t}(t) = £(r) + %Q(t). Then the effective action can be derived as 

S = ( iuj n + H + itsignu n )cl(uj n )c a (uj n ) - Ar](iu n )(n^{-iu n ) - n^-iun)) 

1 



+ E 



+ 



2gA 



2A 2 

-((iuj n )({-iuj n ) + -jj-M(u n + Q )r)(iu n )ri(-iu n ) 



U 



(((iu; n )r}(-iu> n ) + r)(iuj n )((-iuj n )) 



(8) 



where f2 = QJU e s_/U. The effective interaction U e g, is the sum of the Coulomb repulsion 
U and the lattice relaxation energy = g 2 /(2MQ 2 ) as U e s, = U + 4-Elr- The phonon 
Green's function d{iuj n ) = (Q( iuj n )Q{—iuj n )) is given by 

where Xr/(^Yi) = (r](iu n )r](—iu! n )) is the orbital susceptibility of the electronic system. The 
first term is the usual phonon Green's function with the renormalized frequency Q, which 
is higher than the bare Q. This Cl does not depend on the electron response function and 
hence temperature independent. When U <C -Elr ; ^ is much larger than Q and the first 
term of Eq.(9) is irrelevant in the phonon frequency region uj ~ Q ~ O.leV. In the opposite 
limit U ^> -Elr, the renormalization of the phonon frequency is small. The second term 
depends on the temperature (Kohn anomaly) and shows the characteristic lineshape with 
the broadening. From the above consideration together with the experimental fact that no 
temperture dependence of the phonon frequency has been observed, we conclude that the 
Coulomb interaction U is larger than the lattice relaxation energy -Elr- 

Now let us study the electronic system in detail. Following the analysis in section III 



C of ref. ||24j| , the quantum fluctuation between two orbital configurations gives rise to the 
orbital Kondo effect with the Kondo temperature Tk- We obtain an estimate for the Kondo 
temperature as 

RU c s. 



Tk ~ t exp 



A 



(10) 



where 



R = 2x\J 1 + Ay- 1 , for y » l/max(x, x ) (11) 
R = ]jl + jjx 2 (4 + y), for y « l/max(x,x 2 ). (12) 

Here we have introduced the dimensionless variables x = E^/Q and y = U/E^. The main 
conclusion here is that when x > 1, i.e., E^ > fl, the reduction factor R is proportional 
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to x even in the limit U ^> £7lr- This is because the overlap of the phonon wavefunction 
enters the tunneling matrix element even when the Coulomb interaction is dominating. 
This Kondo temperature Tk gives the effective bandwidth B in the large-ci model, and the 



effective mass enhancement is estimated as m*/m = t/Tx ~ exp 



This strong mass 



A 

enhancement manifests itself in the physical quantities as follows. The orbital susceptibility 
Xrii^n ~^ uj + iS) has a peak at uo = with the peak height of the order of t/^/Tk and 
with the width of the order of Tk- The specific heat coefficient 7 should be proportional to 
B~ x , while the coefficient A of T 2 in Eq.(2) scales as A oc B~ 2 . The ferromagnetic transition 
temperature T c is expected to be proportional to B, and hence we expect a large isotope effect 
when R > 1. In other words, the large isotope effect does not rule out the large Coulomb 
interaction U. The edge of the photoemission spectrum at the Fermi energy and the Drude 
weight are reduced by the factor m/m*. As mentioned above, many experiments except 
for those of specific heat seem to suggest that the manganese oxides in their ferromagnetic 
state belong to the strongly correlated region. Then it is natural that the large isotope 
effect is observed. The band width is rather difficult to pin down accurately, but a rough 
estimate from A in eq.(2) is B ~ O.leV, which is comparable to Q. This suggests that 
the suppression is not extremely large. Another experimental fact is that the temperature 



independent Jahn- Teller phonon is observed at the frequency near the original one |L4 
which means that U is larger or at least of the order of -Elr- Then we conclude that the 
relative magnitude of the energy scales is Q < -Elr < U for the ferromagnetic phase in the 
manganese oxides. The small specific heat coefficient necessitates the physical mechanisms 
which have not been included in our model, which will be discussed in the next section. 

Summarizing this section, we have studied a model of manganese oxides in the ferromag- 
netic state taking into account both the Coulomb repulsion U and the Jahn- Teller interaction 
.Elr- These two interactions collaborate to induce the local orbital moments. In the strong 
coupling case, i.e., U e s, = U + 4^lr » t (t : transfer integral), the overlap of the phonon 
wavefunctions (~ e~ ELR ^ n ) enters the tunneling amplitude between the two minima for the 
orbital moment. The phonon spectral function consists of two parts, i.e., the sharp peak 
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at the renormalized frequency (l = ^lyU c s,/U and the broad peak with the width of the 
order of the band width. These results can reconcile the large isotope effect on T c and the 
apparent temperature independent phonon spectrum assuming U > E LR > fi. 

IV. CONCLUSIONS 

We have studied the vital role of the orbitls in manganese oxides. The most impor- 
tant ingredient is the strong Coulombic interaction, which induces the orbital polarization. 
This polarization controls the dimensionality of the electronic system when ordered. In the 
isotropic ferromagnetic state, large orbital fluctuation is expected which gives rise to the 
quantum liquid state below the orbital Kondo temperature Tk- /,From this viewpoint, the 
double exchange interaction is switched on below the orbital Kondo temperature T K , and the 
ferromagnetic transition temperature T c is roughly identified with Tk- This scenario might 
explain the photoemission experiment which shows that the discontinuity at the Fermi level 
grows only below T c . However there still remain several puzzles as follows, (i) the small spe- 
cific heat, (ii) violation of the Kadowaki- Woods relation, (iii) absence of the quasi-particle 
peak in the angle-resolved photoemission spectra in double layered manganese oxides, (iv) 
the residual resistivity larger than the Mott limit. One possible explanation for these is 
the (dynamical) phase separation. The phase separation occurs when the total energy E(x) 
as a function of the hole concentration x has more than two minima and it is possible to 
draw a common tangential line. This is found to occur for several models. Here we want 
to point out that when the bandwidth strongly depends on the orbital configuration, the 
x-dependence of it naturally gives rise to the phase separation. In reality the long range 
Coulomb interaction will suppress the phase separation, but which still remain as a fluc- 
tuation, i.e., dynamical one. Suppose phase separation occurs, and assume that r is the 
metallic fraction while 1 — r is the insulating fraction. Now we relate the observed physical 
quantities to the metallic ones and the fraction r. The T-linear specific heat 7obs^ is reduced 
by the factor r, and hence 7 Q b s = ?"7metai- Concerning the resistivity, the percolating path 
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will lower the capacitance energy when the voltage is applied, and the observed resistivity 
is enhanced by the factor 1/r, i.e., p b s (T) = p me tai(T)/r. Therefore the Kadowaki- Woods 
ratio is related to 



The observed ratio is 10 2 times larger than the universal ratio, and if we attribute this 
discrepancy to the factor r, we obtain r = 0.2. This also suggests 7 mc tai — lAmJ/mol, which 
means that the mass enhancement is ~ 10. Furthermore the Drude wieght is reduced by r in 
addition to the mass enhancement, and is expected to be about 1/50 of the noninteracting 
value in agreement with experiment. The residual resistivity is enhanced by 1/r and can 
be larger than the Mott limit. Roughly speaking, the Hall constant, on the other hand, 
is not modified by the phase separation. The extremely large noise might be also related 
to the dynamical fluctuation of the phase separation. Thus the phase separation seems to 
explain many puzzles in the manganese oxides. However the detailed studies including the 
time/length scale of the fluctuation are left for future investigations. 
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Figure captions 



Figure 1. Mean field phase diagram as a function of the carrier concentration (x) and the 
antiferromagnetic interaction between t 2g spins (Js)- The energy parameters are chosen to 
be a = 8.1 and (3 = 6.67. Schematic orbital structure in the each phase is also shown. 
Dotted line (Js =0.009) well reproduces the change of the spin structure experimentally 
observed. 
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